Hadronic EDMs, the Weinberg Operator, and Light Gluinos 
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We re-examine questions concerning the contribution of the three-gluon Weinberg operator to the 
electric dipole moment of the neutron, and provide several QCD sum rule-based arguments that the 
result is smaller than - but nevertheless consistent with - estimates which invoke naive dimensional 
analysis. We also point out a regime of the MSSM parameter space with light gluinos for which 
this operator provides the dominant contribution to the neutron electric dipole moment due to 
enhancement via the dimension five color electric dipole moment of the gluino. 



I. INTRODUCTION AND SUMMARY 

New sources of CP violation in supersymmetric exten- 
sions of the standard model are highly constrained by 
the null experimental results for the electric dipole mo- 
ments (EDMs) of neutrons and heavy atoms [1,2]. Typ- 
ically, when the superpartners have an electroweak scale 
mass, Avi/, the additional CP violating phases are con- 
strained to be of O(10~^). When confronted with the 
natural expectation that the size of these phases in the 
soft-breaking sector should be of order one, this creates 
a problem for weak-scale supersymmetry (SUSY). 

The interactions which generate EDMs are described 
by a CP-odd effective Lagrangian, induced at IGeV by 
integrating out heavy standard model particles and su- 
perpartners, which contains a series of operators of in- 
creasing dimension. The leading 0-term, 
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has dimension four, and an arbitrary value for 6 consti- 
tutes the usual strong CP problem as its contribution 
to EDMs is unsuppressed by any heavy scale. More- 
over, the existence of additional CP-odd phases in the 
soft-breaking sector of the MSSM aggravates this prob- 
lem by inducing a large additive renormalization of 9 
that survives in the decoupling limit. The conventional 
'cure' - the Peccei-Quinn mechanism - eliminates 9 and 
leaves the dimension five quark EDMs and color EDMs 
(CEDMs), 
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and the Weinberg operator [3] , 
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as the dominant mediators of CP violation from the soft 
breaking sector to the observables. Note that although 
the quark (C)EDMs have dimension five, chiral symme- 
try requires that the corresponding coefficients are pro- 
portional to a light quark mass, and thus di, di, and w 
generically scale in the same way with the overall SUSY 
breaking scale. 

Extracting constraints on the underlying CP-odd 
phases thus requires quantitative knowledge of the depen- 
dence of observable EDMs on , dj , and w normalized at 
the hadronic scale. Recently, the dependence on di and 
di has been determined more precisely using QCD sum 
rules [5], and in this note we turn our attention to the 
Weinberg operator. Although rather intractable within 
the standard framework, we will present several sum- 
rule based estimates. The resulting preferred range for 
the neutron EDM, 



dn{w) = e (10 - 30) MeVw(l GeV) 



(4) 



is a factor of two smaller than conventional estimates 
[3,4] using 'naive dimensional analysis' (NDA) [6]. This 
moderate suppression can be understood through the ap- 
pearance of combinatoric factors which are not accounted 
for within NDA. However, while our result for dn{w) is 
smaller than the NDA estimate, and thus dn{di,di) gen- 
erally dominates the contributions to (i„, there is a regime 
in which d„ (w) is important as it is generated rather dif- 
ferently from the quark (C)EDMs within the MSSM. 

In order to explain this point recall, first of all, that 
there are several generic 'strategies' for curing the SUSY 
CP problem. The first is to require that the superpart- 
ners are heavy enough to suppress all operators of dim> 5 
generated at the SUSY threshold. This decoupling is 
usually applied to sfermions of the first two generations 
only, in order to avoid problems with fine-tuning in the 
Higgs sector. However, this approach is only partially 
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successful as relatively large EDMs may be generated 
through higher loops [7] or through four-ferniion oper- 
ators induced by Higgs exchange [8] . Secondly, one could 
conceive of a universal conspiracy leading to cancellations 
between different contributions [9] , but this is difficult to 
reconcile with the null results for all three types of EDM 
measurement (neutron, paramagnetic and diamagnetic 
atoms) that a priori have different phase dependence [10]. 
A third, perhaps more elegant, option is to invoke an ex- 
act CP or parity at some high-energy scale and specify 
the mechanisms that break supersymmetry in such a way 
that all the relevant soft breaking parameters are ren- 
dered real. This could also be one way of obviating the 
need for axion relaxation [11]. However, some of these 
scenarios may face problems when confronted with the 
large CP violation that is by now well-documented in 
the i?-meson system [12]. 

Given these difficulties, one may pursue another op- 
tion which is to suppress the SUSY contributions by cre- 
ating some (mild) hierarchies between the soft breaking 
parameters in order to suppress the EDMs generated at 
one-loop. Notably, in the limit where gauginos are much 
lighter than the sfermions, all one-loop contributions to 
the EDMs of light quarks and the electron take the fol- 
lowing form: 

di(one loop) ~ (loop factor) x — ^lm(mxA), (5) 

with a similar expression for di induced by the relative 
phase of n and m\. Here i = e,u,d,s, and mgf stands 
for a generic sfermion mass. It is easy to see that as 
mx the expression (5) for di vanishes. Thus a mild 
hierarchy mx ~ (10~^ — 10^^) irisf would appear to be 
sufficient to evade the SUSY CP problem [13,14]. In 
slightly different language, it follows from (5) that in this 
regime the quark EDMs are demoted to dimension seven 
operators and thus are relatively harmless. 

While the quark EDMs are suppressed by this hier- 
archy, we emphasize that sending mx down to hadronic 
scales actually enhances the neutron EDM via the gener- 
ation of the Weinberg operator. The main point is that 
the gluino CEDM, 

jCx = \dx r'^^^a ■ G«75A^ (6) 
can be induced by a top-stop loop [15], leading to 

^2 

dx{one loop) ~ (loop factor) x — |-lm(At — /x* cot /3) 

(7) 

in a basis in which the gluino mass is real. Thus dx 
is a genuine dimension five operator, dx ^ 1/Aw, for 
~ M ~ msf ~ Aw- It follows that for mx ~ Ahadr, 
the gluino takes part in the strong interactions and con- 
tributes to the energy density of all hadrons. Conse- 
quently the neutron EDM is unsuppressed by any addi- 
tional scale, and at a crude level rf„ edx ^ 1/A^y. 
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FIG. 1. Schematic behavior of the neutron EDM dn as a 

function of the gluino mass. Lowering m\ from the SUSY 
threshold there is an initial suppression of dn due to the de- 
crease ofdi{m\) asm\ decreases from Aw to the intermediate 
value m'^* . A further decrease ofm\ in the interval from m^^"^ 
to Ahadr leads to the increase of dn due to the contribution of 
the Weinberg operator, induced by the gluino CEDM. When 
m\ is smaller than Ahadr, dn receives a linear suppression by 
mx. 

This enhancement by the gluino CEDM (6) persists in 
the intermediate hierarchical regime Ahadr ^ mx ^ A^^ 
where, on integrating out the gluino, one generates a con- 
tribution to the Weinberg operator alluded to above that 
scales as l/{mxAw)- At a critical scale mx ~ m™* < A-^ 
these contributions will dominate over di ~ Ahadr^A/A^ 
and dn will start increasing while mx decreases. As we 
will determine below, the scale 

m^"* ~ (6 - 12) GeV (8) 

sets an effective threshold for the maximal suppression of 
EDMs possible with this superpartncr hierarchy*. 

Our results suggest that at this scale the neutron EDM 
is still considerably larger than the experimental bound, 

Ad„(m';^'') ^ (40-80) d^-P, (9) 

unless the SUSY CP phases are fine-tuned. Note that 
both m'"* and rf„(m™*) depend, in addition, on possi- 
ble inter- generational hierarchies for the squark masses. 
When the first generation of sfermions is taken to be 
heavier than A^y, m™* increases while rf„(m™*) de- 
creases. 



*As an aside, we note that (perhaps surprisingly) a gluino 
mass of order m';^"* is still not ruled out by direct constraints, 
and indeed has recently been revived [16] in relation to the 
enhanced hadronic 6-quark production observed at CDF and 
D0 [17]. 
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Therefore, the Weinberg operator has an important 
role to play in minimizing the suppression possible within 
the light gluino regime. Note that for m\ <C Ahadr, the 
CP- violating phase can be rotated to m\ itself leading to 
a suppression of d„ by rriA /Ahadr as one approaches the 
super Yang- Mills limit. A schematic plot of the behaviour 
of dn{mx) is shown in Fig. 1. 

In the next section we turn to the problem of esti- 
mating the contribution to (i„ induced by the Weinberg 
operator, justifying the result (4). Then, in Section 3 we 
describe in more detail the calculation justifying the ar- 
gument outlined above which uses the Weinberg operator 
to limit the suppression of EDMs for light gluinos. 

II. NEUTRON EDM INDUCED BY THE 
WEINBERG OPERATOR 

Unlike the case of d„ induced by the 0-term, or the 
EDMs and CEDMs of quarks, where chiral loop [18] and 
QCD sum rule-based calculations [5] are available, the 
matrix element that relates d„ with the Weinberg oper- 
ator is imknown. The standard estimate, first obtained 
by Weinberg [3], makes use of 'naive dimensional analy- 
sis' [6,19] which keeps track of dimensions, in terms of the 
generic hadronic scale Ahadr, and Goldstone- mediated in- 
teractions through the effective dimensionless coupling 
Ahadr //tt- One finds [3,4], 

dn w{^i) e^QMeY w{ijl), (10) 
47r 

at a low-energy normalization point /i, taking Ahadr ^ 
47r/7r ~ 1.2 GeV. The large value ~ 47r for the coupling 
amounts to demanding that loop corrections are quali- 
tatively similar to the tree-level terms at the matching 
scale. In the gluonic sector, which is important here, 
this means that within the UV quark/gluon descrip- 
tion the relevant value of the gauge coupling is neces- 
sarily very large and consequently the inferred match- 
ing scale doesn't mesh easily with expectations from 
the chiral sector [6,19]. In the present context Wein- 
berg [3], and many papers since [4], have, for the pur- 
pose of evaluating the gauge coupling, chosen a spe- 
cific matching scale corresponding to Qs = 47r/-\/6, or 
as — 2 (cf. as(l GeV) ~ 0.4). If we adopt this normal- 
ization scale in (10), and use (somewhat optimistically) 
the 1-loop anomalous dimension for w [20], the relation 
'^{t^{9s = 47r/\/6)) ~ 0.4 w(l GeV) leads to the most 
commonly used estimate for dn{w): 

di^)~e40MeVu;(/i = lGeV). (11) 

We will avoid quoting a result for the dependence of rf„ on 
w{Kw)^ as there are additional threshold contributions 
from db and dc generated by top-stop-gluino loops, which 
are in general model-dependent [21]. 
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FIG. 2. Perturbative insertion of the Weinberg operator 
into a quark line. The resulting correction to the propagator 
is proportional to w'Y5{qgs{Ga)q) . 

To get some intuition regarding the estimate (11), we 
can consider more carefully the loop-factors which are ef- 
fectively set to unity in (10). For illustration, consider 
reducing the Weinberg operator to the EDM by 'integrat- 
ing out' the gluons. This leads to an effective loop fac- 
tor of gg/{4iT)^ which reproduces (10) provided we take 
Qs ~ 47r. One obtains a similar conclusion for the effec- 
tive scale by considering the gauge kinetic term itself [6] . 
As a consistent matching condition we might then choose 
/U = ^{qs = 47r), leading to a result 

d^f) ~ e 18 MeV u;(m = 1 GeV), (12) 

which is half the size of (11). Although both results (11) 
and (12) are consistent within the expected precision of 
NDA, it is clear that independent quantitative calcula- 
tions are needed to determine dn{w) to better than an 
order of magnitude. 

As a quantitative test of the NDA estimates, we will 
now revisit the calculation of dn{w) using QCD sum 
rules, leading to a result that is a factor of 2 smaller 
than (11) and consistent with (12). To proceed, we note 
first that the leading contribution to the EDM from the 
operator product expansion (OPE) of the nucleon cur- 
rent correlator in the presence of the source (3) exhibits 
a logarithmic infrared divergence. This signals [22] the 
presence of additional operators, required to resolve the 
divergence, whose contributions are generally rather dif- 
ficult to calculate directly. Therefore, we will be content 
to regulate the log-divergent contributions with an IR 
cutoff. These terms will then form the basis of our es- 
timates as they are correspondingly enhanced and thus 
provide the dominant contributions to the EDM. 

We begin by noting that the Weinberg operator allows 
for a perturbative insertion into the quark propagator. 
The leading CP-odd correction is described by the dia- 
gram shown in Fig. 2, and standard manipulations [23] 
lead to the following result: 

iS{p)='-^ + '-^^,{qgs{Gcj)q), (13) 

where the value of the quark-gluon condensate is given 
by [24] 

{qgs{Ga)q) = ml\{qq)\ ~ O.SGeY^ {qq), (14) 
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with (qq) = —(230 McV)^. It is the l/p"^ momentum de- 
pendence in the second term of Eq. (13) which leads to 
the logarithmic infrared divergence ahudcd to above in 
the correlator of two nucleon currents. This signals the 
breakdown of the OPE, but also singles out this insertion 
as providing the dominant effect which we will use in cal- 
culating dn{w). The ambiguity of the infrared log does 
of course render the result less reliable than the corre- 
sponding determination of d„(di,rf,) [5], but nonetheless 
sufficient for our estimates. 

The insertion present in the second term in (13) be- 
haves as a "soft 75-mass". Indeed, while irrelevant for 
large p^, at hadronic scale momenta it mimics the exis- 
tence of an effective CP-odd mass of order 
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[qgs{Ga)q) ~ (120 MeV - 160 MeV)^w 

(15) 



where A^adr is the effective hadronic scale, which we take 
to lie in a range from nip up to ^i^fT^, with previous re- 
sults [5] suggesting that the lower end of this range is 
most relevant for EDM observables. This determines the 
neutron EDM according to the scaling relation [5], 



dn 
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e (1.5-7) MeV «;(Ahadr), (16) 



where the large range in this estimate arises from the 
allowed variation in Ahadr- 

This result is 5-10 times smaller than the conventional 
NDA estimate (11). This is actually not too surprising 
once we recall that a priori dn(w) should be of 0{{qq)^) in 
the chiral limit, while the contribution in (16) is 0{{qq)) 
and thus may indeed be subleading. To test this one can 
consider an explicit sum-rules based estimate [25] utiliz- 
ing the insertion (13). One finds that for the natural 
chirally invariant Lorentz structure, {7^, {Fa)j^} [5], the 
tractable contributions are of 0{{qq)'^) and render a re- 
sult for dn within the range (16). Previous experience 
[5] would suggest that the terms of 0{{qq)'^) are sub- 
dominant, but unfortunately the (a priori) leading con- 
tributions of 0((gg)°) for d„(w) are intractable in this 
direct approach due to the presence of unknown conden- 
sates. 

This analysis suggests that the range (16) might rep- 
resent an underestimate of (i„(w). A natural path to fol- 
low is to consider the sum-rules in chirally-variant chan- 
nels such as (Fa) or ^{Fa)-^ from which one can still 
extract dn{w) along the lines considered previously for 
dn{0) [26]. A convenient means of estimating dn{w) in 
this vein is to calculate the 75-rotation of the nucleon 
wavefunction induced by the Weinberg operator and de- 
termine dn in terms of the corresponding rotation of the 
neutron anomalous magnetic moment /i„: 



dn l^n 



{N\^{GGG)\N) 
mnNi^5N 



(17) 



This approach was considered previously by Bigi and 
Uraltsev [27] who estimated {N\{GGG)\N) in terms of 
{N\GG\N) and the corresponding vacuum condensates. 

We can follow this route and perform a more explicit 
calculation by evaluating the '75' term in the standard 
mass sum-rule correlator of the two nucleon currents. For 
the conventional choice of the loffe interpolating current 
for the neutron 77 [28], we obtain at leading order. 
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d^xe*f^(7?(0)7?(x)) = 
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p'\n{-Aly/p^){qq) 

+ ■■■. (18) 
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It is the relative coefficient between the structures 1 and 
175 that determines the chiral rotation and consequently 
enters into the estimate of [27]. Prom (17) and (18) we 
obtain 

dn ^ Mn^^ti^ln(MVMffl) ^ e 22 MeV w{l GeV), 

(19) 

where we took M/nm = 2 and gs = 2.1. It is impor- 
tant to note that the estimate (19) arises at 0{{qq)'^), 
which we would expect to be dominant, and is indeed 
considerably larger than the estimate (16). A more in- 
volved calculation of the nucleon current correlator in 
an external electromagnetic field [25] reveals additional 
contributions to (i„(u>), but the overall result remains 
quite close to (19). Additional induced corrections, from 
Peccei-Quinn relaxation, would also be subleading [27] as 
they cannot contribute at 0{{qq}^). 

The only other QCD sum-rules estimate of dn{w) that 
we are aware of was made by Khatsimovsky [29] who con- 
sidered a high order term in the OPE proportional to the 
dimension-eight operator F{GGG). An estimate of the 
nonlocal correlator, / d'^x{Q\T{{GGG){Q), {GGG){x)}\Q) 
produced a result for dn{w) similar to (10). However, 
combinatoric factors were ignored in this calculation 
which clearly reduce the result to a value consistent with 
- or somewhat smaller than - (16,19). In practice a pre- 
cise calculation along these lines does not appear feasi- 
ble, as multiple perturbative insertions of the gluon field 
strength into a quark line generally leads to power-like 
infrared divergences [22], signifying the breakdown of the 
OPE. 

Putting these results together, and ignoring the lower 
range of (16) for the reasons discusssed above, we find 
the preferred range for dniw), 



dn{w) ~ e (10 - 30) MeV w(l GeV), 



(20) 



which is a factor of two smaller than the conventional 
NDA estimate (11), and consistent with (12). This re- 
sult will be discussed in more detail elsewhere [25], but 
we turn now to a regime of the SUSY parameter space 
for which this contribution to dn is nonetheless very sig- 
nificant. 
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III. ENHANCEMENT VIA GLUINO COLOR EDM 

As described and schematically illustrated in Section I, 

the neutron EDM is particularly enhanced in the domain 
■'^hadr < rux <C A.w where the gluino develops a color 
EDM via a top-stop loop [15], 



1000 F 



dx{Aw) = - 



tl 




, (21) 



with 5t = Axg[At - /x*cot/3] and fg{y) = (1 - y + 
ln(i/))/(l - yY. Note that this expression is independent 
of m\ and scales as 1/h.w The corresponding contribu- 
tion to the Weinberg operator [3,4,30,31], 



Aw(mA) = — 



"iglimx) dx{mx) 
327r2 mx 



(22) 



scales as l/mx^w- It is worth noting that in addition to 

the obvious enhancement by a factor of Kw /"t,a relative 
to the standard scenario [30], the gluino CEDM-induccd 
shift of the Weinberg operator is also enhanced relative 
to that induced by c or 6 quarks which is of order 1/A^ 
[31,21]. 

The normalization of Aw at the hadronic scale involves 
running the gluino CEDM from down to the gluino 
mass threshold, and subsequent running of w down to 
Ahadr- For completeness, we give the 1-loop /? function 
coefficient, (3q = 11 — 2nx — 2ng/3, where stands 
for the number of light x particles at the scale under 
concern. Besides this, the anomalous dimensions of dx 
and w are given, respectively, by 7"^ = — 18 + /3o and 
= —36 + 3/3o- The latter has been computed in [20], 
and the computation of the former is similar to that of 
the quark color EDMs [32]. 

We now illustrate numerically the impact of light 
gluinos on the neutron EDM using the range for dn{w) 
in (4). Using (21) and (22), we can write 



Mn ~ 100 sin 5t 



(4-12) GeV 



mx 



(23) 



where we have taken M~ = 200 GeV, M~ = 700GeV, 

(h[ = 7i'/4, and the current experimental bound on the 
neutron EDM is d^^'^ < 6 x IQ-^e e.cm [1]. The final 
results are presented in Fig. 3, where we have chosen the 
mid- value c?„ = 20 McV w in (4) . The solid curve stands 
for mx = 1 GeV (with dn/d^^'P « 980 at 5t = 7r/2), 
the dashed curve for mx = mb (with dn/d^^ « 170 at 
6t = 7r/2), and the dot-dashed curve for mx = 20 GeV 
(with dn/dl^'P « 30 at St = tt/2). Thus, for light gluinos, 
where mx ^ (1—4) GeV, one finds that d„(w) exceeds the 
experimental bound by at least two orders of magnitude 
throughout the entire preferred range in (4) unless the 
SUSY phases are tuned such that 5t ^ 10~^. 




FIG. 3. The St dependence of the gluino contribution to 
dn for mx = 1 GeV (solid), mx = mi, (dashed), and 
mx = 20 GeV (dot-dashed). 



Of particular interest is the maximal suppression that 
one can achieve for the EDM in this hierarchical regime 
with light gluinos. We denote by m™* the critical scale 
at which the 1-loop contribution induced by quark EDMs 
(and CEDMs) is approximately equal to the contribution 
associated with the Weinberg operator discussed here. 
Choosing the soft-breaking parameters in the first gen- 
eration of squarks to be 0(200 GeV), and assuming no 
accidental cancellations, we find 



(6 - 12) GeV 



(24) 



accounting for the range in (4), for which the (minimal) 
correction to the EDM is approximately. 



(40- 80)C^ 



(25) 



which still exceeds the experimental bound by at least an 

order of magnitude unless the CP-odd phases are small. 

It is interesting to compare our estimates for dn with 
those one obtains when the gluino is heavy, mx ~ A^. In 
this case, the Weinberg operator is first generated at the 
weak scale at two- loop order [30] . On including the con- 
tributions arising at the 6~quark and c-quark thresholds, 
one finds that dn obtained via (10) only exceeds d'^P by 
at most one order of magnitude [21]. Consequently, the 
light gluino scenario actually induces a larger contribu- 
tion to dn via the color EDM of the gluino. Thus, while it 
is possible to suppress the one-loop contributions to the 
EDMs of leptons and hadrons by taking light gauginos 
[14], the induced contribution to the Weinberg operator 
means that the constraints on the SUSY CP-odd phases 
are not correspondingly relaxed. 

We thank Louis ClavcUi, Glennys Farrar and Oleg 
Lebedev for useful discussions. The work of D.D. was 
supported in part by DOE grant DE-FG02-94ER40823. 
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